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The thermodynamics of the 0{N) nonlinear sigma model in 1+1 dimensions is studied. We 
calculate the pressure to next-to- leading order in the 1/A'^ expansion and show that at this order, 
only the minimum of the effective potential can be rendered finite by temperature-independent 
renormalization. To obtain a finite effective potential away from the minimum requires an arbitrary 
choice of prescription, which implies that the temperature dependence is ambiguous. We show that 
the problem is linked to thermal infrared renormalons. 



I. INTRODUCTION 

The 0{N) nonlinear sigma model (NLSM) in 1+1 
dimensions has been studied extensively at zero temper- 
ature as a toy model for QCD. It is a remarkably rich 
theory, which is asymptotically free and has a dynami- 
cally generated mass gap. It is renormalizable both per- 
turbatively and in the 1/A^ expansion. Moreover, for 
= 3 it has instanton solutions. Unlike the NLSM 
in more than two dimensions, where the theory is no 
longer renormalizable, there is no spontaneous symme- 
try breaking of the global 0(N) symmetry. This reflects 
the Mermin-Wagner-Coleman theorem 1,2], which for- 
bids spontaneous breakdown of a continuous symmetry 
in a homogeneous system in one spatial dimension at 
any temperature. Moreover, the model suffers from in- 
frared (IR) divergences in perturbation theory, since the 
fields are massless in that case |3]. However, it was con- 
jectured by Elitzur Q and shown by David that the 
infrared divergences cancel in 0(iV)-invariant correlation 
functions. In addition, a mass gap is generated nonper- 
turbatively. In the large-A'^ limit, which is equivalent 
to summing all so-called daisy and superdaisy graphs, 
m = /iexp (— 27r/(7^), where g is the coupling constant 
and /X is the renormalization scale. 

Dine and Fischler Q investigated the NLSM in 1 + 1 
dimensions at finite temperature. They calculated the 
free energy in perturbation theory and in the large-iV 
limit. In the weak-coupling expansion, they showed that 
the two-loop contribution to the ideal gas vanishes and 
that the three-loop contribution is infrared finite; the lat- 
ter in fact also vanishes . The leading-order calculation 
in the 1/A^ expansion shows that a thermal mass of or- 
der Ng^T arises. This is a nonperturbative result that 
shows that one is effectively dealing with a gas of massive 
particles. 

In this paper, we extend the analysis of Dine and 
Fischler to next-to-leading-order (NLO) in the 1/A^ ex- 
pansion. At zero temperature, the effective potential (or 
equivalently the Gibbs free energy) has been investigated 



at this order by Biscari et al. "s"]. The l/N correction to 
the thermodynamic potential has the interesting feature 
of containing rcnormalon singularities. We will show that 
it cannot in general be renormalized in a temperature- 
independent way, except at its minimum as a function 
of m? . Away from the minimum, one will need to intro- 
duce a temperature-dependent prescription to deal with 
the poles in the Borel plane. 

Much is known about IR renormalons in the 0{N) 
NLSM in 1+1 dimensions 0, ES IH III El, but the 
consequences at finite temperature have not yet been in- 
vestigated. Thermal renormalons have been studied by 
Loewe and Valenzuela [T^ in cfA theory in 3+1 dimen- 
sions. In this theory, one deals with ultraviolet (UV) 
renormalons only and thus it resembles QED rather than 
QCD. They show that the residues of the UV renormalon 
poles in the Borel plane (which are on the positive real 
axis, whereas in QCD they would be on the negative real 
axis due to asymptotic freedom, such that they do not af- 
fect the Borel transform IS]) in general are temperature 
dependent, but the positions of the poles are not. We 
will show that this is also the case for IR renormalons, 
except at the minimum of the effective potential, where 
also the residues are temperature independent. 

Blaizot et al. ^Bl have recently studied the Gross- 
Neveu model in 1+1 dimensions at finite temperature at 
NLO in the 1/A^ expansion. While there are similarities 
between this model and the NLSM, such as dynamical 
mass generation and asymptotic freedom, no problems 
related to IR renormalons are encountered in Ref. ^3 
(see also E3)- can uniquely define the effective po- 
tential at NLO at nonzero temperature. 

The paper is organized as follows. In Sec. II, we dis- 
cuss the NLSM at zero temperature. In Sec. HI, we calcu- 
late the finite-temperature pressure at NLO. In Sec. IV, 
we discuss various approximations and compare them 
with exact numerical results. In Sec. V, we show that 
one cannot define an off-shell effective potential and this 
is related to thermal infrared renormalons. In Sec. VI, 
we summarize and conclude. 
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II. ZERO TEMPERATURE 

The Euclidean Lagrangian for the nonhnear sigma 
model is 

>C = \{d,'ff + \a{^^-Ng-^), (1) 

where the scalar field $ = (0i, (/)2..., 0jv) forms an A^- 
component vector and a is a Lagrange multiplier that 
enforces the constraint <i>^(a;) — Ng""^. The auxiliary 
field a is now written as the sum of a space-time inde- 
pendent background and a quantum fluctuating field 
a; a — m? + a. The Green's functions of 4> require wave- 
function and coupling constant renormalization in the 

expansion, as discussed by Rim and Weisberger 
through NLO. 

The Lagrangian in Eq. is quadratic in the fields 
$ and the integral over $ can therefore be done exactly. 
One then obtains (cf. e.g. p^ ) 

[ ~ \ N r 2 2-1 
Z = Da exp < ^trln + m + cij + 



N J dr I dx 



^ ,2„-2 I l~„-2 



— m q 
2 ^ 



2«5 



, (2) 



where (3 = 1/T, such that at zero temperature [3 = oo. 
The next step is to expand the functional determinant 
around the classical solution 5 = and integrate over a. 
By scaling a a/y/N, it is seen that this expansion is 
equivalent to a 1/N expansion. 

It is important to realize that is by definition the 
vacuum expectation value of a, i.e. it is the quantity with 
respect to which we will minimize the effective potential 
in order to obtain the pressure. Beyond leading order in 
the 1/N expansion receives divergent contributions. 
In order to show this, we rewrite the expression for m? 
in terms of m|, where is defined as the pole of the 
propagator D^{P,m) of <!>. At NLO, D^{P,m) is given 

by 



where ^ = ^1 + Am^/P^. Choosing renormalization 
point P^ = —mi., one obtains |2nl |. 



(6) 



where li(x) is the logarithmic integral, 

\i{x)^V j\t^^. (7) 

Here, A denotes the ultraviolet momentum cutoff and V 
indicates a principal-value prescription for the integral. 
Solving Eq. © for m? , we obtain 



2 2 / A2 

m = m i — h — 5" 



(8) 



In order to have a pole with a residue equal to unity, one 
needs = 1 — jjJ^'{P, m)\p2^_^2 , which yields 



N 



(9) 



The wavefunction renormalization Eq. 10 is in accor- 
dance with that obtained by Flyvbjerg ,20j. Rim and 
Weisberger flj calculated the wavefunction renormal- 
ization constant in dimensional regularization and it 
also agrees with Eq. © as can be seen by identifying 
ln{A'^ /fi'^) 2/e where d=2-e. 

The effective potential V through next-to-leading or- 
der in the 1 /N expansion is given by 



V = 



- -N 
2 

d^P 



d^p 

(2^ 



In 



In [n(P, m) 



(10) 



where we have added a subscript h to g to indicate ex- 
plicitly that it is the bare coupling constant. Evaluating 
the integrals in Eq. H1U|) using an ultraviolet momentum 
cutoff A, one obtains 



D^{P,m) 



P^ + m^- jf^{P,m) 



(3) 



where Z^p is the wavefunction renormalization constant 
and 



E(P, m) 



d^Q 



1 



(27r)2 (P + Q)2+to2 n(Q,m) 



(4) 



is the self-energy function. Here, H(P, m) is the inverse 
propagator for a: 



I f d^Q 1 



1 



)2 Q2 + ^2 (p _^ QY + ^2 



(5) 



V 



'^9t 
1 

-TO h I — - 



8tt 



1 -f In- 



A2 



TO^ 



(A2 



A2 

2to2) In In — - 

1 2 



TO/ 



2to2 



7B 



1 - In 



A2 \ 
4m2 / 



(11) 



Here and in the subsequent results we have dropped to- 
independent divergences and terms that vanish in the 
limit A2 ^ CX3. 

To obtain the pressure, one evaluates the effective 
potential at its minimum. The condition for the mini- 
mum is given by equation 



dV 
dim? 







(12) 
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Eq. H12() is often referred to as a gap equation. Differen- 
tiating Eq. 1)11(1 . one obtains 



47r 



2\ A2 1 / A2 
1 1 In — H I 21nln — 



N m2 N 



A2 

-li— + 27£; + 41n2 



(13) 



To see that Eq. 1(13(1 becomes finite after coupling con- 
stant renormalization, we first express it in terms of m^, 
using Eq. ©: 



An ( 2 \ A2 2 / 

^7 = 1 In — 7T H In In — ^ 

al \ N) ml N\ ml 

+"fE -f- In 4 



(14) 



The renormahzation constant for g is denoted by Z^^ 
and is given by 



r I 2 \ A2 1 g2 A 2 

1 + i- I 1 - — I hi — — ^Inln — 
47r 



NJ fi^ N2tt 



(15) 



Making the substitution g'^ Zg2g'^(^), we obtain the 
renormahzed gap equation: 



The expression Eq. ((T31) for Z^^ is exact in g^{fi) up to 
order l/iV^ corrections and results in the known NLO 



/3-function |all8|,|21|: 

Pid) - aM = _(i_^)£^(1 + _£^) ,(17) 



dA 

'^^^ ^ ^ d/^ 27r V N 



N J 271 
4 / 2 



7V27r 



(18) 



Using the gap equation, one can obtain the value of 
the effective potential V at the minimum, where it equals 
the pressure V. In terms of bare quantities, we obtain 



1 



= -(iV-2) A^ln 



An 



(19) 



This equation will be used to subtract the pressure at zero 
temperature from the pressure at finite temperature. 



III. FINITE TEMPERATURE 

The results at zero temperature are obtained ana- 
lytically, but at finite temperature this is in general not 
possible. Therefore, we will investigate the pressure nu- 
merically. However, we are able to isolate the ultraviolet 
divergences analytically. 



The effective potential through next-to-leading order 
in 1/N is now given by 



m'^N 



In [n^(F,m)] 



(20) 



where the inverse propagator 11^ (P, m) is the finite tem- 
perature version of Eq. |(SJ) and we have defined the sum- 
integral 



= T 



E 



dp 
2^ 



(21) 



Summing over Matsubara frequencies and averaging over 
angles, 11-^ (P,m) reduces to 

n^(P,m) = r ^R{P,q)[l + 2n{E,)] , (22) 

where Eq — \/ q"^ + m^ and n{x) = (exp(/3a;) — 1)^^ is 
the Bose-Einstcin distribution. The function P(P, q) is 
given by 



R{P,q) 



1 



P2 + 2pq 



An (P2 + 2pqY + AplE^ 



(23) 



The inverse propagator cannot be evaluated analytically, 
so we will evaluate it numerically. For this purpose, it is 
necessary to isolate ultraviolet divergences analytically. 
As expected on general grounds, i.e. from the absence 
of temperature-dependent ultraviolet divergences, and as 
verified numerically, the quantity 



Pi = ^ Inn^ (P,m)-y i^lnn^(P,m), 



(24) 



is finite. To calculate Pi we used an Abel-Plana for- 
mula p3 | . In order to isolate the divergences we consider 
the limit p 3> T, where we can approximate 



lf{P,m) « n(P,TO) 



1 



p2 



An P4 + Am^pl '^^ 



= nSE(P,m) 



where 11 (P, m) is given in Eq. and 

dq 



JiiPm) 



E„ 



n(Eq 



(25) 



(26) 



In order to split off the prefactor of the 
logarithm in 11^^*', we define n(P, m) = 
-47r A/p2(p2 + 4m2)n(P, m). This gives the following 
contribution to the free energy 



d^P 



In [p2 



{2nY 

2 / ^2 



Am 



= 1 + In , , 

2n \ Am^ J 



(27) 
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In order to isolate the infinities, we need the large-P be- 
havior of {P,m): 



The calculation of the self-energy E(P, m) at pq 
and = —mi, at finite temperature yields 



p2 



2m^ 



^k^{P.m) = In— — (1 + Ji) 



(28) 



where 



m^exp(— Ji). This yields 



where 



d^P 
(2^ 



■In 



ngE(P, "^) = £'2 + finite terms , (29) 



1 

47r 



Inln- 



A2 



m li ■ 



A2 



m 
'2^ 



In In ■ 



A^ 



Finally, we define 

F2 - 



d^P 
(2^ 



■In 



n(p, m) 



Do 



(30) 



(31) 



Again we have checked numerically that the quantity P2 
is finite, demonstrating that we have identified all ultra- 
violet divergences. 

Putting everything together, the finite temperature 
effective potential becomes 



V 



Nm" Nm- 



1 + ln- 



A2 



N , 
8^ 



(Pi + D1+F2 + D2) 



where 



8 /"^ 



dqq^ 

~e7 



(32) 



(33) 



We again note that we have systematically dropped m- 
indcpendent divergences and terms that vanish in the 
limit A^ — > cxD. 

The gap equation H12() at nonzero temperature now 
becomes 



477 A2 1 

~9 = 111 + 



A2 

2 Inln — 



-2 In 



A^ 
4m^ 



- 47r 



TO^ dm^ 
d(Pi + P2) 



dm^ , A^ 



dm^ 



(34) 



From the fact that is temperature independent, one 
can conclude that is also temperature independent at 
leading order in the expansion, when it is a solu- 
tion to the gap equation. We will use this fact later on 
to conclude that the pressure can be renormalized in a 
temperature-independent way. 



ml - 



N 



,. A2 

mi 



(35) 



where P3 is a finite function that depends on the temper- 
ature as well as m<^. Since we use m^ merely as a way 
to express the renormalized gap equation in terms of fi- 
nite quantities, any choice of P3 will do and we choose 
P3 = 0. We have checked numerically that other choices 
indeed do not alter the final result for the pressure. Using 
Eq. Eq. (PH) now becomes 



47r 



= In 



A2 



1 

N 



2 In In 



A2 



-21n- 



4m? 



- 47r 



d(Pi + P2) 
dml 



(36) 



To render the gap equation finite, we again only need 



to make the substitution 



Zg^g ifJ-), where Z 2 is 



9b ' '^g^f \t-/: 
given by Eq. (|15|1 . The renormalized gap equation then 
becomes 



47r 



2 

N 



N 



Ji(/3to0) +ln4 + 27r 



d(Pi + P2) 
dml 



(37) 



Using the gap equation, we obtain the value of the 
effective potential at the minimum which is equal to 
the pressure. Using Eq. H35|) and expanding the Jq and 
Ji functions, one ultimately obtains for the pressure at 
nonzero temperature minus the pressure Eq. (|19|l at zero 
temperature. 



V = — V 



T=0 



N -2 



(0) - ml] 



N 
8^ 
1 
2 



+ ^[T^MPm4,)+mlJi{pm^)] 



2 d(Pi+P2) 
mj, : — 7t i<i — i<2 



dmi 



where Pi and P2 are functions of and m? 



m' 



and ?Ti^(0) = m^iT = 0). We have numerically evaluated 
the expression for the pressure, after solving Eq. (|37|l for 
m^{T). The result for different values of N is shown in 
Fig. n for the arbitrary choice g^{fi = 500) = 10, hence 
T is given in the same units as /z. As can be shown and 
seen in the figure, V/NT'^ approaches an A^-dependent 
constant (to be evaluated below) at large temperatures, 
which for iV — *■ 00 is 7r/6. Moreover, it approaches zero 
in the limit of zero temperature. If we normalize the 
pressure, for a given value of N, to its value at T = 00, 
we find that the normalized pressure has a very small 
dependence on N. 
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FIG. 1: Pressure as a function of temperature at NLO for 
different values of A*'. 



constant is 



z-^ = 1 

9^ 



El 

in 



1 - 



N 



In- 



A2 



(42) 



which is consistent with the perturbative renormahzation 
constant to leading order in g^. Making the substitution 



7^ Zg2g'^, we obtain 



47riV 



TV Ji - Kt - + 2{N - 2) In ■ 



where the function is 



Kf = ±4 



Ap uj± 
mEj, 



n{uj±) 



(43) 



(44) 



IV. HIGH-TEMPERATURE 
APPROXIMATIONS 

In Ref. [23, Bochkarev and Kapusta consider the 
nonlinear sigma model in 3+1 dimensions, which is non- 
rcnormalizable, at NLO in the 1/N expansion. Since the 
result for the pressure cannot be obtained analytically, 
they resort to a "high-energy approximation". We will 
make the same approximation and compare it with the 
exact numerical results obtained in Sec. III. 

The idea of the high-energy approximation is that in 
the part of 11-^ proportional to the distribution function 
n{Eq), the important contribution comes from the region 
where po,p 3> q. One can therefore approximate the 
self-energy Ji^{P,m) by its high-energy behavior. In the 
present case, this amounts to 



n'^(F,m) w n(P,m) 



1 



Ji 



(39) 



where u!± = y^p'^ + m? ± m. This expression is identical 
to Eq. H25|) . After simply discarding the T = contri- 
bution to n^(P, m), as done in Ref. (l^, the effective 
potential is approximately given by 

+^5? In W + "+] + 5^ In [pS + . (40) 

The resulting expression for the gap equation is 

i 1 



Ng- 



P2 + m2 4f 

-i 



P niEp pI + Lo\ 



mEp pI+uj^ 



(41) 



Again, the gap equation requires coupling constant renor- 
mahzation. In this approximation, the renormahzation 



Note, however, that the pressure is finite even when 
we substitute the unrenormalized gap equation H41|l into 
Eq. gn|): 



1 

8^ 



[{K+ +K^)T^ + {K+ -2)n?] , (45) 



where the function is 



y2 



dp p^ 

En 



n{uj±) 



(46) 
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FIG. 2: Pressure as a function of temperature at NLO for 
different values of A'^ compared with the high-energy and high- 
temperature approximations. 

From Fig. [3 one can see that the high-energy ap- 
proximation underestimates the pressure compared to 
the exact result. The advantage of an approximation like 
the high-energy approximation is that the analytic cal- 
culations are simpler and that it is easier to implement 
numerically. 

We suggest a different approximation, which is bet- 
ter than the high-energy approximation. We will calcu- 
late the inverse a-propagator 11"^ by first integrating over 



6 



the momentum. We obtain 



1 



(47) 



In the limit m <^ T, we can approximate H^{P,m) by 
n^rp(P, m), where keep only the qn = mode in the sum 
Eq. g7I): 



1 1 



2 l3mP^ + Am'^p^ 



(48) 



Since it follows from the leading order gap equation that 
for high temperature and for all values of the coupling 
constant, m <C T, we call this approximation the high- 
temperature (HT) approximation. By using that + 
4m^p^ can be written as [pq + {p + imf + + (p ~ 

jm)^ + vn?] and shifting p — > p ± im after taking the 
logarithm, the functions Fi and F2 can be approximated 
by 



ZTT 







(49) 



A numerical calculation of i^i + F2 shows that for 
m/T this approximation has an error smaller than 
10 percent. Approximating Fi and F2 using the high- 
energy approximation is less accurate. The result for the 
pressure in the high-temperature approximation is shown 
for comparison in Fig. |3 (again for (7^(/Lt = 500) = 10). 

One can approximate the pressure even further by 
expanding the functions Jq and Ji in the limit (im —> 0: 



Jo 



Ji 



A-Kfim - 2 log 



4^ 

+0{{Pmf) , 
— + 2 log^^ 



(3m 
47r 



- lE 



(50) 

0{{Pmf) . (51) 



Inserting the approximations given in Eqs. (|49f) and l|51f) 
into the gap equation H37|) . one obtains 



(3m 



2tt 

7m 



In 4 



1 + TT 



7b 



in 



(52) 

which indicates that (3m ~ 1/lnT for large T. In the 
limit m/T — > 0, we obtain for the high-temperature ap- 
proximation of the pressure 



m 

1^ > (53) 



7VT2 
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V N 



where the first term is the pressure of a gas of free mass- 
less particles with — 1 degrees of freedom. 



V. THERMAL RENORMALONS 

We have shown that a finite pressure at finite tem- 
perature can be obtained after subtraction of the zero- 
temperature pressure and coupling constant renormal- 
ization. This agrees with the general expectation that 
ultraviolet divergences are connected with short-distance 
physics and therefore independent of the temperature. 
While we have shown this explicitly at NLO in the 1/A^ 
expansion, this is not the case for the effective potential 
away from its minimum. 

In the expression Eq. Hll|l for the effective po- 
tential at zero temperature, the two contributions 
A^ lnln(A^/rn^) and m^li(A^/m^) cannot be removed us- 
ing m- independent counterterms ^. While this may not 
be a problem at zero temperature, it would certainly 
become one at finite temperature when m becomes a 
function of temperature. This would imply temperature- 
dependent renormalization, which is not acceptable. In 
Ref. 0, these two divergences are dealt with by con- 
sidering the effective potential normalized to its zero- 
mass value, i.e. Vnlo(™) — Vnlo(O). This subtraction 
is ill-defined due to infrared divergences and therefore 
one should understand it as subtracting the contribu- 
tions from ln[n(P, m)] obtained in the limit P^/m^ 
00 24] . This is called the "perturbative tail" . If we 
denote Il(P,m) in this limit by IloQ{P,m), one finds 
U^{P,m) = ln(PVm2)/(47rP2) and 



d^P 

^ln[n^(P,m)] 



1 

An 



A^ In In 



A2 



-TO li 



A^ 

TO^ 



(54) 



where we have implicitly used the principal-value pre- 
scription. In Refs. 0, IS ll^, this subtraction is not 
motivated, but we point out that the subtracted con- 
tribution is associated with IR renormalons. As shown 
in jljj, the vacuum expectation value of a, i.e. to^, is in- 
herently ambiguous, when one tries to separate (in order 
to subtract) perturbative contributions proportional to 
A^ from the nonperturbative ones proportional to in 
the limit where A^ ^ m^. In 11], it was shown that 



(OjajO) 



-I- 



47rTOg^o 
N 



d^p 1 an 
J2^iid^ 



(55) 

where to^q = A^ exp(— 47r/gj) and the 1/A^ contribution 
arises from the tadpole diagram shown in Fig. |31 

One can show that this equation is in agreement with 



the gap equation H13|) if we write ' 



TO 



LO ■ 



' TO 



NLO 



/N. 



The part of the integral in Eq. H55|) that has the IR 



1 Note that the quantity In In — In In ^ (with fj. ^ m) 
diverges as A'^ ^ 00, whereas In In A^. _ In In ^ vanishes. 
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FIG. 3: Tadpole diagram contributing to (a) at next-to- 
leading order in 1/A'^. The wavy line represents the a propa- 
gator and the solid line the $ propagator. 



renormalon pole in the Borel plane is in fact the contri- 
bution from the integrand in the limit P^/m^ —^ oo: 

^ d?P 1 OToo _ ^ U 

(27r)2 ^ ^4^ 

1 A2 

= -1 2e""Ei(a;), (56) 

47r m2 

where x = \iy[A? / m^) . In the limit x —> oo, the logarith- 
mic integral has the asymptotic expansion: 



ambiguity, we have also not considered obtaining a fi- 
nite effective potential or even a finite pressure at zero 
temperature. However, we have calculated the quantity 
V"^ — 7^"^^°, which is free of renormalon ambiguities and 
is finite after temperature-independent coupling constant 
renormalization 

Finally, we comment on the possible temperature 
dependence of renormalon contributions to (a) and the 
effective potential. One can show that Eq. at finite 
temperature has exactly the same renormalon contribu- 
tion, i.e. neither the pole nor the residue become temper- 
ature dependent. Secondly, the perturbative tail of the 
effective potential which is given by the first two terms 
of D2, corresponds to poles in the Borel plane at 6 = 
and 6=1. The positions of the renormalon poles are 
not affected by temperature. Only the residues become 
temperature dependent, except at the minimum of the 
potential, as we concluded earlier. The fact that renor- 
malon pole positions are not affected by temperature, but 
residues are, is also the case for the thermal ultraviolet 
renormalons in in 3 1 dimensions studied by the 
authors of Ref . ,1^ . 



e-^Ei(a;) 



=p tire 



n=0 
poo 



db 
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(57) 



where arg(6) — ±e. From Eq. H57|l . it is clear that there 
is a renormalon pole at 6 = 1. This shows that when A —> 
00 the value of (0|a|0) is inherently ambiguous at NLO, 
due to the freedom in the choice of prescription. David 
has shown that this ambiguity also arises in dimensional 
regularization ,.1QJ . 

The same problem appears in the calculation of the 
effective potential, but not in the gap equation. The 
latter can be seen from the last term of Eq. H1U|I which 
contributes to the gap equation as follows 

1^ 1 /-^^ 

2 9m2 J (27r)2 2 J {2n)^ U Om^ ' ^ ' 

The ambiguity that would arise from this term when re- 
moving its perturbative tail (cf. Eq. H56|l ) cancels in the 
gap equation l(T^ against the one arising in (cf. Eq. 

The perturbative tail of the effective potential, i.e. 
the first two terms of D2 defined in Eq. (|30|l . corresponds 
to poles in the Borel plane at 6 = and 6=1, respec- 
tively. Since fh? is only temperature independent at the 
minimum (at LO only, but that is sufficient since we are 
working at NLO), the subtraction of the perturbative 
tail will become temperature dependent, except at the 
minimum. Since subtracting temperature-dependent di- 
vergences renders the remaining temperature-dependent 
terms ambiguous, we refrain from following this strategy 
and thus from trying to define a finite effective potential 
at finite temperature. In order to avoid any renormalon 



VI. SUMMARY AND CONCLUSIONS 

To summarize, we have calculated the pressure in 
the NLSM at finite temperature to NLO in the ex- 
pansion. Our main result is that we obtain an unambigu- 
ous, finite pressure, by subtracting the zero-temperature 
value of the pressure and renormalization of the coupling 
constant in a temperature-independent way. This proce- 
dure cannot be carried out away from the minimum of 
the effective potential and we have argued that defining 
a finite, effective potential by the subtraction of the so- 
called perturbative tail, leads to ambiguities associated 
with IR renormalons. In general, these become temper- 
ature dependent, and this casts doubt on the usefulness 
of defining a finite effective potential. 

We have calculated the expression for the pressure 
at finite temperature numerically and observe that the 
1 /N expansion is a meaningful expansion for all temper- 
atures. We have also investigated the high-energy ap- 
proximation that was originally applied to the NLSM in 
3-1-1 dimensions by Bochkarev and Kapusta. In 1-1-1 di- 
mension, where one can compare with exact numerical re- 
sults, we have shown that it underestimates the pressure 
for all temperatures. We have suggested an improved 
approximation, the so-called high-temperature approxi- 
mation. This approximation has the advantage that it is 
quite easy to produce numerical results and agrees better 
with the exact results. At asymptotically high tempera- 
tures the pressure approaches that of a gas of iV — 1 free 
massless particles. 
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